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Abstract. We calculate exactly the time-dependent reduced density matrix 
for the central spin in the central-spin model with homogeneous Heisenberg 
couplings. Therefrom, the dynamics and the entanglement entropy of the central 
spin are obtained. A rich variety of behaviors is found, depending on the initial 
state of the bath spins. For an initially unpolarized uncntanglcd bath, the 
polarization of the central spin decays to zero in the thermodynamic limit, while 
its entanglement entropy becomes maximal. On the other hand, if the unpolarized 
environment is initially in an cigenstate of the total bath spin, the central spin 
and the entanglement entropy exhibit persistent monochromatic large-amplitude 
oscillations. This raises the question to what extent entanglement of the bath 
spins prevents decoherence of the central spin. 



1. Introduction 

The central-spin model, or Gaudin model, describes one spin coupled to A — 1 bath 
spins via isotropic Heisenberg interactions, 
w-i 



including a magnetic field which couples to the central spin only. Our focus here is on 
spin- 1/2 particles. 

The importance of this model is manifold. From a quantum information 
processing point of view, it describes the interaction of one qubit with an environment, 
serving as a model to study decoherence processes. In this regard, it is used to capture 
the hyperfine interaction of an electron trapped in a quantum dot with the nuclear 
spins [1, 2, 3, 4, 5, 6, 7, 8, 9]. In [3, 4], it was found that non-uniformities in the 
couplings Aj induce a non-exponential decay of the central spin. This result was 
obtained by a perturbative calculation for zero polarization of the bath spins, and by 
an exact approach in the thermodynamic limit for a fully polarized bath. 

Apart from that, from a theoretical point of view, the challenging feature of this 
model is the following: The model is exactly solvable by Bethe ansatz (BA)[10]. On 
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the other hand, in the thermodynamic limit, exact results can be found in a mean- 
field approach [11, 9], where quantum mechanical operators are replaced by their 
expectation values. Especially the compatibility of the quantum-mechanical solution 
for a large number of particles with the mean-field solution was addressed very recently 
[9]. There the dynamical evolution of the central spin in the quantum and mean- 
field solutions was observed to be similar for zero magnetic field and homogeneous 
couplings. That observation was made with the bath being initially in a spin-coherent 
state with a fixed total spin. The work summarized in the present article shows 
explicitly that such a similarity between the classical and quantum solutions is not 
possible if the bath is initially in a product state without entanglement between the 
bath spins. 

Despite the exact solvability of the model, far away from the thermodynamic 
limit, no systematic approach for dealing with the entire spectrum of the Gaudin 
model is known. Richardson [12] gave a method to calculate finite-size corrections for 
the pairing Hamiltonian, whose BA equations are equivalent to those of the Gaudin 
model. Explicit results for finite-size corrections to the ground state and to low- 
lying excited states for the pairing Hamiltonian have been obtained [13] starting from 
Richardson's work. However, these techniques have not been applied to the spectrum 
of the Gaudin model. The only resort seems to be numerics [1, 2, 14, 15], which is 
hard enough a problem in itself. 

Given the importance of the Gaudin model from both theoretical and applied 
points of view it is highly desirable to obtain analytical results for an arbitrary 
number of particles N with an arbitrary number of flipped spins M (as compared 
to the ferromagnetic "all up" state) which corresponds to a total magnetization 
S z = N/2 — M. In this work, we consider the special case of homogeneous couplings 
and zero external field, Aj = 2Vj and h = (in [4], this model was considered for 
M = 1, 2). For this choice of parameters, the time-dependent reduced density matrix 
Po(t) of the central spin is calculated. Therefrom, the quantum mechanical average 
(<So}(f) = tr {Sgpo(t)} and the entanglement entropy E(t) = — tr {po{t) log 2 Po(t)} 
are derived. In all cases considered here, the central spin and the bath are uncorrelated 
initially, with the central spin in the "down" -state | and the bath state specified 
further below. 

Although it seems as if real physical systems like quantum dots are not described 
by homogeneous couplings [2], the study of that simple special case turns out to be 
very fruitful. First of all it is solvable explicitly for arbitrary N, M, which enables us 
to present exact results for finite particle number and arbitrary polarization. Indeed, 
highly non-trivial results, such as the decoherence time or the amount of screening 
of the central spin, could be obtained and are presented below. Furthermore, it is 
important to study the influence of different factors on the dynamics separately in 
order to assess their respective importance. These factors include the limit of large 
particle number, the magnetization and entanglement of the bath, the non-uniformity 
of the couplings and the magnetic field. The focus in this work is on the first three 
of these factors, with the couplings homogeneous and without a field. Finally, the 
study presented here is expected to be helpful for the investigation of generalizations 
(namely, to non-uniform couplings and finite h- field). Exact results for these more 
general settings are highly desirable from an applications point of view, and the study 
of the homogeneous case will serve as a point of reference. 

Our main results are summarized in the following (units are chosen such that 
h = 1). The simple structure of the Hamiltonian leads to an eigenvalue spectrum in 
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which all transition frequencies (energy differences) are commensurate. Therefore all 
physical quantities are strictly periodic with period Tp, the Poincare recurrence time. 
The recurrence time does not depend on the system size; Tp — tt (2tt) for N even 
(odd). The detailed dynamics within every period Tp depends on details of the initial 
state of the bath spins; the initial condition for the central spin is such that it always 
starts in the "down" state. 

We first discuss the case of an unpolarizcd (or minimally polarized) and 
unentangled bath state; that is, the bath contains M — 1 « N/2 — 1 down spins 
and N — M up spins in a product state. We find that the Fourier spectra of both 
the spin expectation value (S§)(t) and the entanglement entropy E(t) contain M 
frequencies. However, these frequencies combine such that (S§)(t) approaches zero 
within a decoherence time t c i ~ 1/y/N. Accordingly, E(t) reaches its maximum 
value 1 within a time 0(l/y/~N). Deviations from those limits are 0(1/N). In other 
words, in this case, the bath provides a perfect screening of the central spin in the 
thermodynamic limit. 

Next we still consider a bath state with zero or minimal polarization, but with 
entanglement among the bath spins, namely an eigenstate of the total bath spin 
operator. In contrast to the previous case, now both (Sq)(£) and E(t) oscillate with 
only one frequency and maximum amplitude. The frequency depends on the total 
bath spin quantum number. Since the initial eigenstate of the bath is a multipartite 
entangled state, this raises the question to what extent an initial bath entanglement 
generally protects the central spin from decohering. We do not answer this question 
in the present article, but rather leave it as an interesting perspective for future work. 

We turn to a partially or fully polarized bath. Then in the time evolution from 
an unentangled initial bath state again M frequencies are present, which combine 
such that the dominating ones are 0(N — 2M), that is, at the lower end of the 
frequency spectrum. In the thermodynamic limit, (Sg)(t) = —1/2 + 0(1/N), 
E(t) = 0((\nN)/N), and even/odd effects are 0(1/N). The leading-order non- 
constant terms oscillate increasingly rapidly as N grows, see Fig. 3 below. 

Finally, the partially polarized bath can be initially in a total spin eigenstate such 
that the dynamics of (S§)(t) and E(t) are again driven by one frequency only, around 
a mean value and with an amplitude which depend on the total bath spin quantum 
number. 

When comparing these results with the behaviour obtained from a mean-field 
approach, we find similarity only with the evolution from an initial eigenstate of the 
bath. This is consistent with [9] in so far as there, similarities between the quantum 
and mean-field approaches were observed for an initially entangled bath. However, 
the evolution from the unpolarizcd initially unentangled product state has no classical 
analogue. This result is worthwile being stated: Although the exact classical solution 
of the model (1) is known [11], the question which initial bath states reproduce this 
solution on a quantum-mechanical level has not been addressed yet. 

This article is organized as follows. In the next section, we briefly sketch our 
approach for calculating p (t) and present the dynamical evolution of (S§)(t) and 
E(t). We distinguish between an initially unentangled environment, the environment 
being initially in an eigenstate and the mean-field solution. Technical details of the 
calculation are deferred to an appendix. The article ends with an outlook on possible 
further research directions. 
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2. Results 



For Aj = 2Vj and h = 0, the Hamiltonian can be written in terms of conserved 
quantities as follows: 

H = S 2 -S 2 b -Sl (2) 

where the bath spin Sb = &j commutes with the Hamiltonian as does Sq. The 

total spin is defined as S = Sb + Sq. The Hamiltonian has eigenvalues 

F _ [ S b , S = S b + 1/2 

1 -S h -1, S = S b - 1/2 ' 

with S b = (N- 1)/2,(N- l)/2 - 1,. .., N/2 - AI + 1/2. 

Our aim is to calculate (S(j(t)) for a fixed value of the total magnetization 
S z = N/2 — M (in this case, the same symbol is used for the operator and its 
eigenvalue), starting with an initial configuration | JJ.)| . . .)(,. We consider the two 
cases where the initial bath state | . . .)& is a pure product state (not entangled) and 
where the bath is entangled initially. Our approach consists in decomposing the initial 
state |-||)| . . .}b into a sum over energy eigenstates \<p) with constant magnetization S z , 

\il)\...) b = J2cM 

<t> 

where <j> stands for the quantum numbers S b ,S — Sb ± 1/2. Starting from this 
decomposition, the density matrix p(t) is written down and therefrom the reduced 
density matrix po(t) := trb {p(t)} of the central spin is calculated, where the bath 
degree of freedoms are traced out. For any further details on the calculation of po(t) 
the reader is referred to appendix Appendix A. 



2.1. Initial product state 

For an initial product state | JJ-; |, ■ • • , |, 1\ • • • 1 1") with N particles, Mb down spins in 
the bath, M = Mb + 1 down spins in total, 2M < N, the reduced density matrix for 
the central spin is calculated in appendix Appendix A, Eqs. (A. 8, A. 9). The results 
for po(t) and (S z )(t) are rewritten here in a slightly different way: 

p ° (t) = (N~u)l^-k)\L {2{M ~ k)(N ~ M ~ m cos [{N ~ 2k)m){t 1 



k=0 



+ ((M - kf + (N - M - kf + 2(M - k)(N - M - k) cos [(TV - 2k)t}) |J|)(J||} (3) 

who- ^e>-*w-*» 



2 ^ (N — 2k){N — k)\k\ 

x |(JV — 2M) 2 + 4(iV - M - k)(M - k) cos [{N - 2k)t]} . (4) 

The entanglement entropy is given by E(t) = — tr {po(t) log 2 po{t)}. The case 
2M > N is related to the case 2M < N by the symmetry property 

(S$)(t)\ M = {S5)(t)\*-M+i (2M>N), (5) 
in other words, the expectation value of the central spin is invariant under a reversal 
of all bath spins at t = 0, with the exception of the case M = N, where the system 
is in a stationary state. The evolution from the initial state | ft; j, . . . , |, j, . . . , |) 
(N particles, M up spins in total) is obtained from Eq. (3) by exchanging ft and JJ-, 
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resulting in an overall sign in Eq. (4). Let us now discuss the results (3, 4) in more 
detail. 

Generally the recurrence time is Tp = ir for even N and Tp = 2tt for odd iV.J For 
a fully polarized bath, i.e. M = 1, only one frequency iV is present, and the amplitude 
of oscillations is 0(1/N). This case as well as the case M = 2 have been considered 
in [4]. In the other extreme, for even N and magnetization S z — N/2 — M = 0, the 
series (4) can be summed for t = ir/2, yielding 

(S z ) (n/2)\ N=2M = 2{N \ i} ~ ± N - oo (6) 

N N 

E{n/2)\ N=2M - log 2 (iV - 1) - 2(iv _ 1) log 2 - 
AT-2 /TV \ 

iog 2 hr - i (7) 



2(N-1) oz \2 
1 AT- 



21n(2)7V 2; 

Applying Stirling's formula to the factorials in (4) for N — > oo, M = pN/2, < p < 1, 
shows that the amplitudes of all but the low-frequency Fourier components are 
suppressed exponentially. Especially, in the unmagnetized case, the sum in (4) can 
be evaluated by the saddle point method for large particle numbers and short times, 
i C rp. To leading order this yields 

/oo 
\x\e~ 2x2 cos(2xViV<)dx, t < r P . (8) 
-oo 

The approximation leading to (8) is commonly called the quasi-static approximation 
(QSA). Thus in the limit N 3> 1 and for times which are short compared to the 
recurrence time, the particle number only enters through the time scale. We will see 
below that this time scale plays the role of a decoherence time. 

We emphasize again that Eq. (8) describes the time-evolution starting with a pure 
product state of the whole system. Previous studies [16, 17] concentrated on an initial 
configuration which involves an average over all bath states and performed the whole 
calculation in the QSA. Our result (4), however, is more general: It is valid without 
any restrictions, neither on the particle number nor on the time, and describes the 
time evolution from a single product state. As a special case, we may choose the 
mixed state of [16, 17] as initial state by carrying out the corresponding average. In 
this case, we recover the results of [16, 17] in the QSA limit. 

Further insight is gained by plotting (S§)(t), E(t) as functions of t for different 
values of the bath magnetization. In Fig. la) the case of even N = 2M is depicted. 
The function (5g)(t) in the interval < t < n is symmetric with respect to t = ir/2, 
showing two zeros and two local maxima. The maxima of (iSoH^) are found to behave 
as {Sg)(t max ) = 0.14 + 0.33//V with t max = 2.1/Vn,tt - 2.1/VN. The zeros scale 
like t = 1.3/VN,7r - 1.3/ y/N. Consequently, (S$)(T d < t < it - r d ) -> for 
= 2M — ► oo, with ~ l/\/~N. Thus Td plays the role of a decoherence time, 
as anticipated above. 

The entanglement entropy reflects this behaviour: Any deviation of |(>S , o)(t)| from 
its maximal value 1/2 results in a finite entanglement entropy of the central spin [1]. 



\ "generally" means that there are special choices of iV, M such that rp is shorter than tt. For 
example, for M = 1, one has Tp = 2-kv/N , v integer. 
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Figure 1. (color online) Central spin expectation value (5q)(0 < t < tt) in an 
unpolarized or minimally polarized bath for JV = 2M in a) and JV = 2M + 1 in 
b), M = 5, 10, 20, 50, 100, 200 starting with an unentangled state of the bath. The 
dots at t = n/2 in a) are the result (6). The inset c) shows the Fourier spectrum 
for M = 200, N = 400 (red circles) and JV = 401 (blue dots), as a function of the 
frequency JV - 2k, k = 0, . . . , M - 1. 



Especially, a zero in (S^)(t) corresponds to a point where E(t) takes its maximal 
value 1, and a local extremum in (Sq)(*) translates into a local minimum for E(t), see 
Fig. 2a). This means that the scaling of the zeros and maxima of (5g)(i) induces an 
analogous scaling of the local extrema of E(t), according to Eqs. (A. 10, A. 11). 

Similar features are found for the case of odd N = 2M + 1, Fig. lb). In this 
case, however, (S§(t)) reaches an overall maximum at t = tt, with (S§(tt)) — > 1/2 as 
N = 2AI + 1 — > oo. The width of this maximum (estimated by the zeros next to it) 
scales ~ 1/y/N, which means that the time interval affected by even/odd effects (in N) 
is O (1/y/N). The maximum (Sq)(tt) translates into a local minimum E(n) — ► with 
N = 2M+1 — * oo, cf. Fig. 2b). To conclude our analysis for the case of an unpolarized 
bath, the Fourier spectrum, i.e. the coefficients of the oscillating terms in (<Sg(t)), 
Eq. (4), as a function of N — 2k are plotted in Fig. lc). Only low-lying frequencies 
contribute, high frequencies arc suppressed exponentially as observed above. 

We now consider a partially polarized bath, M = pN/2. Again by a numerical 
finite-size scaling for a variety of choices of p, we find that the amplitude of the 
oscillation of (S§)(t) at t = n/2 is 0(1/ N), around a mean value -1/2 + 0(1/ N), 
independently of p. Thus for a partially polarized bath, (S§)(t) = (Sq)(0) + 0(1/ N) 
as N -> oo. In agreement with this, E(t) = 0((hxN)/N). Figs. 3a), 3b) show (Sfi)(t), 
E(t) for the special case N = 4M, or p = 1/2. 

Note that all these results are valid for homogeneous couplings. The quantum- 
mechanical time-evolution of (5o)(i) for inhomogeneous couplings has been addressed 
in [3, 4]. In that case, oscillations with frequencies 0(N), (0(y/N)) and 
amplitudes 0(1/N), (0(1)) set in for the initially unentangled polarized (unentangled 
unpolarized) environment, before the broadening of the couplings leads to a non- 
exponential decay. That decay at long times does not depend on the particle number. 
There is no such long-time decay for homogeneous couplings. Instead we observe 
periodic behaviour with a finite recurrence time rp independent of the particle 
number. For a polarized or partially polarized, initially unentangled bath state, the 
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Figure 2. Entanglement entropy E(0 < t < 7r) in an unpolarized or minimally 
polarized bath for N = 2M in a) and TV = 2M+1 in b), M = 5, 10, 20, 50, 100, 200. 
The dots in a) at t = tt/2 are the result (7). In all cases, the bath is initially 
unentangled. 




Figure 3. (color online) Time evolution of a) the central spin expectation 
value (S z (t)) and b) the entanglement entropy E(t) for partial bath polarization, 
N = 4M, M = 5, 12, 25, 50, respectively. In all cases, the bath state is initially 
unentangled. 



oscillation frequencies and amplitudes are of the same order of magnitude as they are 
for inhomogencous couplings [3, 4] during the initial phase before the long-time decay. 
However, the time evolution starting with an initially unpolarized, unentangled bath 
as discussed here is governed by a decay after a short time ~ 1/ \/N, as opposed 
to the long-time decay found in [3, 4]. Interestingly, Td here scales with the particle 
number in the same way as the period of the initial oscillation in [3, 4] docs. 



2.2. Initial bath eigenstate 

Let the bath initially be in an eigenstate of S% (with eigenvalue Sb(Sb + 1)) and with 
magnetization S§ = N/2 — M + 1/2. In this state, the bath is entangled initially 
except for full polarization. 
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The reduced density matrix of the central spin for this initially entangled bath 
state is given in Eqs. (A. 5, A. 6) in appendix Appendix A. From these expressions, 
the time evolution of (Sq) is deduced 

+4(5 fc - S z b + l)(S b + S z b ) cos [(2S b + l)t]} . (9) 

The most striking difference to (4) is that only one frequency is involved in Eq. (9). 
Let us consider again the cases of a polarized, partially polarized and unpolarized 
bath. For a polarized bath one has S§ = Sb = {N — l)/2. In this case, Eqs. (4,9) 
yield identical results, which is to be expected because, as mentioned above, a fully 
polarized bath is in a product state. 

For a partially magnetized bath, S b = pN/2 with < p < 1. This docs not 
fix Sb which can take values S b < Sb < (N — l)/2. For S b = S b , we obtain 
(S§)(t) = —1/2 + 0(1/N), which is similar to the time evolution starting with an 
initially unentanglcd state. It is also this behaviour that one might expect from an 
intuitive point of view: a macroscopic polarization of the bath constitutes an effective 
gap O(N) seen by the central spin, which prevents decoherence. 

However, if we consider the opposite limiting case when the bath is initially in an 
eigenstate with Sb = (N — l)/2, then, in the leading order in N, 

{£%}(£)= - \ {f + (1 -p 2 ) cos Nt), 

which is an oscillation around the mean value —p 2 /2 with the frequency N, a behaviour 
qualitatively different from the one at full bath polarization, Sb = S b . 

Finally, the difference between an initial product state and an initial eigenstate 
of the bath becomes most pronounced in the case of small magnetization 0(1). Above 
(see for example Fig. 1) we found that an initial product state of the bath causes the 
central spin to decohere after a time t<i ~ 1/y/N. However, starting with an initial 
bath eigenstate, the time evolution reads (S z )(t) = — 1/2 cos [(2Sb + l)t] for even N 
and S§ = 1/2. That is, 

{-l/2cos2i, St = 1/2 

: : , (10) 

-1/2 cos Nt, S b = (N-l)/2 
whereas for odd N and S§ = 0, (S§)(t) = - 2(2S l +1) 2 (1 + ^S b {S b + 1) cos [(2S b + l)t]): 

( -1/2, S 5 = 

( - MW) , S b = (N-l)/2 

Thus (5g)(t) oscillates at maximum amplitude with one of (=j) possible 

frequencies (depending on Sb) for odd (even) N. For odd N and S b = 0, the limiting 
bath eigenstates are the completely symmetric (maximal S b ) and antisymmetric 
(minimal Sb) states. Note that in the latter case, the complete system is in a zero- 
energy eigenstate, which means that the time evolution is trivial. Especially, for 
N = 3, these are the S z = triplet and singlet states, respectively, two of the two- 
qubit Bell states. 

The behaviour of the entanglement entropy E(t) is found in a straightforward 
manner from the reduced density matrix in Eqs. (A. 5, A. 6). The qualitative behaviour 
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can be read off from that of (S^)(t) by remembering that (Sq) = ±1/2 corresponds to 
E = 0, (Sq) = is equivalent to E = 1 and any local extrema of (Sq)(£) are translated 
into local minima of E(t). 

2.3. Classical/ Mean- field solution 

The Gaudin model can be treated in a mean-field approach in the thermodynamic 
limit [11, 9], However, as shown here, when the dynamical evolution of the central 
spin is considered, the equivalence between the quantum-mechanical and the classical 
picture does not hold if the initial state is a product state. It does hold, though, if 
the bath is initially in an cigcnstate as is shown in the following. 

The mean-field approach consists in replacing the effective field seen by each spin 
with its quantum-mechanical average [11, 9]. We denote the corresponding quantities 
by an additional subscript MF, 

_^ j 2So • Sb,MF 
\ 2Sq,MF ■ Sb 

Since the resulting equations of motion So.fa = i [H, Sq^] = 2<So,i> x Sb.a.MF are linear 
in the quantum mechanical operators, expectation values are taken here as well and 
one arrives at the classical equations of motion 

So,MF = 2So,MF X Sb.MF 
Sb.MF = 2Sb.MF X So,MF- 

From these equations it follows that 2(Sq.mf + Sb,MF) '■= 9 = const., as expected. 
Let us choose the coordinate system such that g = (0, g, 0) T , and the initial conditions 
S 0i mf(0) = (0,0,S ) T , S b ,MF{0) = (0, 5 /2,-5 ) T . Then 

(— Sq sin gt \ 

S cos gt J 

Comparing this result with the quantum-mechanical solution starting from an initial 
bath eigenstate (10), one notices that both results are closely related, where the 
parameter g/2 in the classical solution plays the role of Sb in the quantum- mechanical 
picture. Note that in both cases, the polarization of the bath along the z-axis is of 
the same order as the central spin. However, for an initially uncntangled bath state, 
the quantum-mechanical solution does not match the mean-field behaviour. 

This leads to the question to what extent initial entanglement in the bath is 
necessary in order to reproduce the classical solution within the quantum-mechanical 
framework. In [9] it was observed that if a spin-coherent (i.e. entangled) state is 
chosen as the initial bath state, the quantum-mechanical solution is recovered by the 
mean-field solution. Given that the exact classical solution of the model (1) has been 
obtained recently [11], addressing that question quantitatively remains an interesting 
task for the future. 

3. Conclusion and outlook 

We have studied the time evolution of the central spin in the central-spin model with 
homogeneous couplings and without magnetic field by calculating the corresponding 
reduced density matrix. Whereas the mean-field solution is an appropriate description 



Spin- dynamics in the central spin model 



10 



for the evolution starting with the bath being unpolarizcd and initially in an cigcnstate, 
the full quantum mechanical solution is needed for an initially uncntanglcd unpolarizcd 
bath. In this case, the central spin is completely screened by the bath spins after 
a decoherence time Td ~ 1/y/N until n — Td, within deviations ~ 1/N for a non- 
polarized bath in the thermodynamic limit. The central-spin polarization recurs after 
the recurrence time rp = ir(2ir) for N even (odd). A complete polarization reversal 
at rp/2 = 7r occurs for N odd; however, the duration of this reversal is 2rd, so 
that in this respect, even/odd-effects are O (1/y/N). The dynamics of the entan 
glement entropy of the central spin reflects the dynamics of the central spin: For 
an unpolarizcd initial product state, the central spin finds itself maximally entangled 
with the bath spins after a time O (1/y/N). Thus for an initially unentangled bath, 
non-trivial dynamics are obtained only from the full quantum-mechanical solution at 
finite particle number - in the thermodynamic limit the decoherence time vanishes, 
and after that, the dynamics is frozen. 

Directions for future research include the investigation of the effect of initial 
entanglement within the bath. Two questions are of particular importance in this 
context: (i) Whether it is generally true that initial bath entanglement protects the 
central spin from decoherence (the eigenstates of a non-polarized bath considered in 
section II. B are entangled - but at this stage it is unclear whether it is rather the 
eigenstate- or the entanglement-property that is important). A hint to a positive 
answer to this question can be found in [18], where it was argued that persistent 
entanglement in the bath serves as a barrier against decoherence of the central spin. 
The second question (ii) is under which conditions initial bath entanglement allows to 
recover the mean-field solution. Furthermore, the exact generalization of the quantum 
mechanical solution to inhomogeneous couplings is highly desirable. First results have 
been obtained [19] in the form of an exact formula for (So)(t) stemming from the 
Bcthc Ansatz. The explicit evaluation of this formula especially for a non-polarized 
bath is the subject of future work. 
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Appendix A. Calculation of the reduced density matrix 

We first have to decompose the initial state 

|4;i,...,U,...,T) == im---4,t,---,T>iv-i 

with N particles and M = Mb + 1 down-spins, 2M < N, into eigenstates of H. 

For a given value of the bath spin Sb, there are two values of the total spin, 
namely S = Sb + 1/2 and S = Sb — 1/2. We denote an A r -particle eigenstate of S 2 
(with eigenvalue S(S + 1)), S 2 (with eigenvalue Sb(Sb + 1)) and S z (with eigenvalue 
S z ) as \S, Sb,S z ) n. Furthermore, an A^-particle eigenstate of S 2 and S z is denoted by 
\S, S z )n- The total angular momentum eigenstates of an arbitrary angular momentum 
plus a spin 1/2 are a standard example (cf. for example [20]) in quantum mechanics 
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texts. For the present case one obtains that [20] 



Sb 



Sb 



, Sb, S b 



— , Sb, Sb 



N 



N 



V2S~b 



y/S b + SZ\1t)\Sb,SS-l)lf- 



+Vs~ b 



V2Sb 



s z b 



l\$)\Sb,S z b ) N - 



y/S h -Sl + l\ti)\S h ,Sl-l) N - 



VSb 



\Sb, Si 



b/N-1 



(A.l) 



(A.2) 



are eigenstates of S b with eigenvalue Sb(Sb + 1) and of H as well with eigenvalues Sb, 
— Sb — 1, respectively. Note that the states on the right-hand side of these equations are 
product states between the bath (which is specified by its total spin and the total spin's 
z-component) and the central spin. The eigenstates (A.l, A. 2) were also constructed 
in [21], where the partition function of the model was studied. The relations (A.l, A. 2) 
can be inverted to obtain 

1 



m\s b ,st) N _ 



m\s b ,si-i) 



^2S b 



1 



^s b - s§ + 1 



Si, 



y/Sb + S% 



Sb 2 ' ^b, Sb 



N-1 



VWb 



\fS~b 



Vs b -s§ + i 



Si, 



1 



1 

h 2 

2/jv. 
Sb, S b 



, Sb, S b 



^)n 



N 



S b —, Sb, S b 



N 



(A.3) 



(A.4) 



If the initial state of the bath is an eigenstate of S b and S b (and hence an entangled 
state in all but the fully polarized cases), (A.3) can be used to construct the time- 
dependent reduced density matrix po(t) for the central spin. The projection operator 
|<f>(t)) (<f>(t)| developing from the initial state of the full system contains combinations 
of the Hamiltonian eigenstates (A.l, A. 2). From that projection operator, po(t) is 
obtained as po(t) = tr& {|$(t))($(t)|}. Thus the trace over the bath spins for each of 
those combinations has to be performed; for example (A.l) yields 



tr. 



Si, 



1 



- , Sb, Sf 



N 



Si, 



-, Sb, S b 



N 



i(Sb + s z b )\t)(l\\ + (Sb-s z b + m)(H 



2S b + l 

and similar expressions for the other combinations. The final result is 
p (f)=a(i)|W| + (l-a(i))|W| 



(A.5) 



a(t) 



2(S b 



SI + l)(S b + SI) 



(2S b + l) 2 

x (1 - cos[(25 6 + l)t}), (A.6) 

implying that for an initial eigenstate of the total bath spin all physical quantities 
oscillate harmonically with frequency 2^ + 1 . In order to find p (0 f° r the case of an 
initially uncntangled bath state, we proceed as follows: Starting from the completely 
polarized state 1 1, . . . , |)at_m, the decomposition of | |l) 1 1, . . . , |)jv-m can be found 
as above. We then consider this state as the new bath state | J., f, . . . , f) n-m+i and 
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find the decomposition of | Jj.) 1 1, j, . . . , T)n-a/+i- This procedure is repeated until the 
state I, . . . , I, |, . . . , |)iv-i is reached. The result is of the following form 




„(N,M) 



AI 



N 



(A.7) 



with unknown real coefficients c 



(N,M) 



The symbol denotes the binomial coefficient 
N\/{M\(N — M)\). Note that the total number of states is 2 • £)£ =( 



i qM-i _ 



as expected. All states are mutually orthogonal, so that in the expressions for the 

2 



r~tlVL - 

.7 = 1 



enters and it is this 



density matrix, the sum of squares dk,N,M 
quantity that we set up a recursion relation for. The time-dependent reduced density 
matrix po(t) corresponding to the state (A.7) is 

Po(f) = /?(f)IWK(i-/3(*)MW 



from which we derive 

(S z (t)) 
E{t) 



dk. 



N.AI 



M 

; g(iV-2fcF 



(M - k)(N - M - k) 



x (1 - cos[(7V - 2k)t}), 



(A., 



(A.9) 



(A.10) 



= -p(t)]og 2 p(t) 

-(l-/3W)log 2 [l-/3(i)]. (A.ll) 

Hence, for an initial product state of the bath (with M — 1 reversed bath spins) a 
spectrum of M + 1 frequencies N — 2k (k = 0, . . . , M) governs the dynamics. In order 
to further determine the function (3{t) (A.9), we next decompose the (A^ — l)-spin bath 
part of the state (A.7) into a single down spin and a (A^ — 2)-spin state with (M — 1) 
down spins in a way completely analogous to (A.7), with the replacements N — > N — 1, 
M — > M — 1, and calculate | -||)| I, ■ ■ ■ , I, f > • • • , T)iV-i again. By comparing coefficients 
in (Sg(i)), the following recursion relation is found for dk,N,M' 

M-l-k 



dk 



N.AI 



N - 1 - 2k 
N-M- 



dk,N-l,M-l 
- (jfe-1 



dk-i,N- 



1,M-1 



with di 



O.N.M 



[c N M -\Y 

dk,N,M 



N - 1 -2(fc- 1) 
This recursion relation is solved by 
(M - 1)\(N - M)\(N - 2k) 



(N - k)lk\ 

which leads to the result Eq. (3). Note that X^fclo 1 dk,N,M = lj which is the 
normalization condition for the state Eq. (A.7) and at the same time guarantees 
(S »)(0) = -l/2. 
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The calculation for 2M > N proceeds analogously. Now the state | 4)1 Tj • • • j T 
! I) • • • i I)jv-i w ith N — M f-spins is built successively starting from the polarized 
state | |, . . . , 1)m-i- I n eacn step, one "[-spin is added according to Eq. (A. 4), until 
finally. | -1J.) is included using Eq. (A. 3). It tuns out that in this case the dynamics of 
the central spin is the same as in the case of N + 1 — M down spins and M — 1 up 
spins in the bath. In other words, (S§}(t) is invariant under reversal of all bath spins 
except for the case of a fully polarized bath, where reversal of all bath spins leads to 
a stationary state of the entire system. This leads to Eq. (5) in the main text. 
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